Radial fall of a test particle onto an evaporating black hole 
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£NJ ■ Abstract 

O I' A test particle falling into a classical black hole crosses the event horizon and ends up in the 

singularity within finite cigcntime. In the 'more realistic' case of a 'classical' evaporating black hole, 
an observer falling onto a black hole observes a sudden evaporation of the hole. This illustrates 
the fact that the discussion of the classical process commonly found in the literature may become 
obsolete when the black hole has a finite lifetime. The situation is basically the same for more 
complex cases, e.g. where a particle collides with two merging black holes. It should be pointed 
out that the model used in this paper is mainly of academic interest, since the description of the 
physics near a black hole horizon still presents a difficult problem which is not yet fully understood, 
but our model provides a valuable possibility for students to enter the interesting field of black hole 
physics and to perform numerical calculations of their own which are not very involved from the 
■ computational point of view. 

O ■ PACS numbers: 04.25.-g; 04.70.-s; 04.70.Dy 

The Schwarzschild metric generated by an uncharged non-rotating classical black hole is given by 

U 

ds 2 = (i_ r JLy dt 2 _ (i _ Hy'dr 2 - rW, (1) 

•rH , 

X ' where r s = 2GM /c 2 is the Schwarzschild radius, G the gravitational constant, Mq the mass of the 
black hole, c the speed of light, and dQ 2 denotes the line element squared on a two-dimensional unit 
sphere. Due to the Hawking radiation, the black hole evaporates, if we exclude any accretion of 
energy. For macroscopic black holes, the luminosity is proportional to the inverse mass squared. The 
corresponding ansatz for the time-dependent black hole mass M{t) 

M{t) = -kM(t)' 2 , (2) 

leads directly to 

M(tf = |(to-t), (3) 

where to = 3Mq//c is the lifetime of the black hole. Accordingly, we have 

r s (t) = k'(t -t) 1 / 3 . (4) 

It is clear that limited knowledge is available about the actual evaporation process of black holes 
at present, which is based mainly on the classic works of Unruh and Hawking PJE1> and there is an 
ongoing interest in the subject [1103113. Therefore, our model is purely academic, but we consider it 
worthwhile to investigate it due to its simplicity and because it exhibits the interesting feature that a 
particle falling onto a black hole does not vanish from our universe, since it experiences an immediate 
black hole evaporation. A realistic treatment of the process presented here would involve detailed 
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knowledge about the physics in the very close vicinity of the black hole horizon. E.g., Babichev et al. 
[H] discussed the case of black holes accreting dark energy in the form of a phantom energy, where the 
energy density p > and the negative pressure p < fulfill the condition pc 2 +p < 0. For a black hole 
surrounded by a constant phantom energy bath, the black hole mass decreases like (t + t) -1 , where 
r > is a characteristic evolution time. Therefore, we could have chosen decay laws for our discussion 
given below that differ from eq. <jl|). but the qualitative results would nevertheless remain the same. 
It is also clear that equation Q will fail to describe the evaporation process of the black hole in its 
final stage. 

We also point out that an introduction to the description of evaporating black holes can be found 
in a known paper by William Hiscock jjj. There, different decreasing functions M(t) in a slightly 
different framework using the so-called Vaidya metric jSJ |§] were considered. The main drawback of 
using simple approaches like those in [3| is the fact that the impact of the emitted radiation is not 
taken into account, because the black hole does not evaporate due to the emission of radiation but by 
artificially adding negative energy to the black hole. 

We give here a handwaving argument in order to illustrate the strange physical conditions that 
might be expected in the vicinity of the black hole horizon. For a distant observer, a massive black 
hole is a very cold object with Hawking temperature 

where kp is the Boltzmann constant. It is well-known that the condition for thermal equilibrium in 
a static system is y/goo(r)T(r) =const., where T(r) is the temperature measured by a local static 
observer. This condition can be most easily derived when two systems in thermal equilibrium are 
considered, which are coupled only via their thermal radiation, such that the black body radiation of 
the two systems has to travel through the gravitational field. For the region near the horizon with 
r = r s + 5, where 5 > is small, goo is given approximately by 



Furthermore, the spatial distance d to the horizon is 




(7) 



Therefore, a hypothetical static observer very close to the horizon with d = lp = ^HG/c 3 measures a 
temperature of the order of the Planck temperature Tp, since 

k B T BH he 3 2r s Tic 

— = k B Tp. (8) 



goo(6 = l p /4r s ) l P l P 

For an external observer, the test particle gets arbitrarily close to the horizon within finite time, but 
in contrast to our considerations made above, it is not in a static state, but accelerated, and does 
therefore not experience the extreme Unruh effect as a static observer. Nevertheless, it will be difficult 
to make any reasonable predictions about the fate of a particle approaching the horizon, where the 
physics becomes Planckian. 

In the following, we assume that the metric generated by the evaporating black hole is given by 
expression but with r s replaced by the time-dependent expression given by eq. Q. In view of 
the very slow evaporation process and the low energy density of the emitted Hawking radiation, this 
is certainly a satisfactory approximation for our purpose. Near the black hole horizon, the validity of 
the approximation is dubious at best. 
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From now on, we will use geometric units (G=c=l) throughout. If we parametrize the (geodesic) 
worldline of a massive test particle by its eigentime r, then the equation of motion of the particle 
falling radially onto the black hole can be derived from the Lagrangian 



C = B(r,t)i 2 -£(r,i)~V, B(r,t) = I 1 



1^2 



(9) 



where the dot denotes differentiation with respect to the eigentime r. The Euler-Lagrange equations 

6C 5C 



Or 



5(i,f) 6(r,t) 



then read 



Eq. (jnjl leads to 



d T (2B(r, t)t) = d t B i 2 - d t {B- x ) f 2 , 
d T (2B{r,t)~ 1 r) = -d r B i 2 + d r (£ -1 ) f 2 . 



d r B . 9 d T B . ■ B n _ . g 

-!— f 2 + — ft d r Bt 2 , 

2B B 2 ' 



, d r B t B 2 dtjB- 1 ) . 2 

t = =-rt —t — — r . 

B 2B 2B 

Additionally, we can take into account that the four-velocity of the test particle is constant 

Bi 2 - B- X f 2 = 1, 



(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



since we used the eigentime r as curve parameter. 

We solved eq. (|12|) in conjunction with eq. (|15|) numerically. We considered first the analytically 
solvable case of a particle falling radially on a stable black hole (with mass M) starting at rest at a 
distance R 

dv 

= 0, r(r = 0) = R. (16) 

r=R 



dT 



Introducing the parameter r/ defined by 



(I- 1 )- 



7] = arccos ( 



(17) 



one obtains 



R 



1 + cos 77), r 



i? 3 



t = 2M\ 



JL 

2M 



1 



R 

V + J^(V + sin 7?) 



8M 
+ 2Mln 



(77 + sin 77), 



g 

2M 



1 + tan ■ 



_R 1 _ t aT1 2 

2M L< * 2 



(18) 



Such a particle reaches the horizon after an eigentime r H 

"It 




R 



arccos -f - 1 \+2\ - -5r 



(19) 



This value served as a practical independent check for the accuracy of our numerical calculation. 

A simple explicit finite difference Euler forward scheme is fully sufficient for the numerical inte- 
gration of eq. H13|) . where one can eliminate t by using eq. Q15|). But it is necessary to adapt the 
eigentime step At during the integration process in order to obtain an accurate solution. As long as 
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Figure 1: Test particle falling into a stable black hole with Schwarzschild radius r s = 1. This graph 
can be found in many standard textbooks. 

the particle has a relatively large distance to the horizon, the time step can be held constant, but in 
the vicinity of the horizon, a update of the time step according to 

Ar n+1 = —Ar n (20) 

r n 

for each integration step is a good choice. The effect is a reduction of the time step when the radial 
coordinate r starts to change faster with respect to the eigentime r, and the resolution of the trajectory 
with respect to the external observer time t is maintained. 

Fig. 1 shows the trajectory of the particle falling into a stable black hole with r s = 1 and R = 2 
with respect to the particle's eigentime r and the external observer time t. In this case, the particle 
reaches the horizon for r H = V2(ir/2 + 1) ~ 3.6357 and the singularity for r s = ^^R 6 jr s ~ 4.4429. 

Fig. 2 shows the typical situation for an unstable black hole. In order to obtain reasonable plots, 
we used 'short' black hole lifetimes, but the qualitative picture remains the same independently from 
the actual lifetime of the black hole which may be very large for massive holes. From the external 
observer's point of view, the particle impinges on the black hole horizon, where it remains captured 
until the black hole evaporates. In the final evaporation phase, the particle even looses contact to the 
horizon due to the fast decrease of the black hole radius. 

From the point of view of the test particle, the black hole seems to evaporate as soon as the particle 
approaches the horizon. The amusing point in our toy model is the fact that from the 'test particle's 
point of view', the evaporation time of the black hole is larger for black holes with a short lifetime. 
This is due to the fact that the evaporation process continuously reduces the distance between the 
horizon to the particle, such that the particle experiences a smaller time contraction near the horizon, 
when the horizon shrinks faster. For massive black holes, the black hole evaporation time approaches 
the time that the particle would need to reach the horizon of a stable black hole with r s = r s (r = 0), 
given by eq. Fig. 3 shows T(io) — t h , i.e. basically the eigentime of the particle with the same 

initial conditions as in Fig. 2 at which the evaporation process of the black hole with initial lifetime 
to comes to the end. 

We conclude with the remark that within our simplified framework one should expect basically the 
same result for more complex processes, e.g. for a test particle approaching two merging black holes 
which evaporate subsequently. The test particle cannot cross the classical horizon, and an asymptotic 
observer will never observe that the particle impinges the horizon within finite time |10| . 




Figure 2: Typical trajectories for a test particle falling onto a black hole with initial Schwarzschild 
radius r s (0) = 1 and a lifetime to=20. 




Figure 3: Eigentime r(io) reduced by the eigentime r H needed by a test particle to reach the horizon 
of a stable black hole, versus the black hole lifetime to. 
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